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We explore the possibility that a new interaction between muons and protons is responsible for 
the discrepancy between the CODATA value of the proton radius and the value deduced from 
the measurement of the Lamb shift in muonic hydrogen. We show that a new force carrier with 
roughly MeV mass can account for the observed energy shift as well as the discrepancy in the muon 
anomalous magnetic moment. However, measurements in other systems constrain the couplings 
to electrons and neutrons to be suppressed relative to the couplings to muons and protons, which 
seems challenging from a theoretical point of view. One can nevertheless make predictions for energy 
shifts in muonic deuterium, muonic helium, and true muonium under the assumption that the new 
particle couples dominantly to muons and protons. 

PACS numbers: 



I. INTRODUCTION 

A recent publication [T] announced a measurement of 
the Lamb shift in muonic hydrogen that seems to require 
a value of the proton's radius, r p — 0.84184(67) fm, which 
differs by five standard deviations from the value given 
in the CODATA compilation 0, r p = 0.8768(69) fm. 
It is of course possible that the reason for the different 
values of r p extracted from muonic and electronic hy- 
drogen lies within the standard model, involving subtle 
QED and/or hadronic effects. It is nevertheless worth 
considering whether new physics could be the explana- 
tion. Here we explore the possibility that the discrep- 
ancy arises from the existence of a roughly MeV-mass 
force carrier that couples the muon to the proton. As we 
will see, such a force can also resolve the long-standing 
discrepancy between theory and observation in (g — 2)^ 
measurements. However, atomic precision measurements 
and neutron scattering experiments at low energies con- 
strain the coupling of this new force to electrons and neu- 
trons, respectively. It is not difficult to construct models 
that effectively decouple the force carrier from the elec- 
tron (by arranging for the particle to couple to mass) or 
the neutron (by arranging for the particle to couple to 
charge), but decoupling it from both is a more serious 
challenge to model building attempts. Taking the force 
to act on muons and protons only, it is still possible to 
make predictions for energy shifts in related systems - 
muonic deuterium, muonic helium, and true muonium - 
with minimal model-dependence. 



II. CONTRIBUTIONS TO THE LAMB SHIFT 
IN HYDROGENIC SYSTEMS 

Independent of whether the new force is mediated by 
a scalar or a vector-boson, the non-relativistic potential 
between the proton and muon can be written as 



where m$ is the mass of the force carrier, s is its spin, and 
<7 M and g p are its couplings to the muon and proton, re- 
spectively. A similar expression holds for the electronic- 
hydrogen system. Note, however, that the sign of the 
potential may be different depending on the relative sign 
between g e and g^. 

This potential gives an additional contribution to the 
Lamb shift in the 2S1/2-2P3/2 transition, which using 
first-order perturbation theory is given by 
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with f(x) = x 4 /(l + x) A . Here — (am MP ) _1 is the 
Bohr radius of the system to leading order and m MP is 
the reduced mass of the /i-p system. A similar expression 
holds for the e-p system, but the Bohr radius is a factor ps 
m^/rrie larger. This expression is convenient for a direct 
comparison with the leading order contribution from the 
proton radius, 
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where r p is the proton's radius, and n(= 2) is the princi- 
ple quantum number. 

If everything else is the same between the e-p and [i- 
p systems, then since f(x) is a monotonically increasing 
function that asymptotes to unity at large values of x, the 
resulting energy shift in the e-p system is always larger 
than the corresponding shift in the /i-p system [3] . If the 
force is attractive in both systems, the apparent proton 
radius will always appear smaller in the e-p system, con- 
trary to observations. An attractive force must therefore 
couple more strongly to muons than to electrons if it is 
to explain the discrepancy in the proton-radius determi- 
nation. 

A different possibility is that g e is not suppressed rel- 
ative to g^, but the force is repulsive, leading to a larger 
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apparent proton-radius in the e-p system. This possibil- 
ity is consistent with the (g — 2) e constraint discussed 
in the next section. However, since the effects of such 
a light force are suppressed at higher momentum trans- 
fer, this possibility seems in tension with the value of the 
proton radius extracted from scattering data which gen- 
erally also imply a larger proton radius 0J. Therefore, 
for the purpose of this paper, we concentrate on the pos- 
sibility that it is a new attractive force that modifies the 
jj,-p system. 

III. CONTRIBUTIONS TO (g 2) e , M 

The scalar and vector-boson contributions to the elec- 
tron and muon anomalous magnetic moments are 16] , 

Aa ' = f (-) a £("W™<), (4) 
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where mi is the mass of the electron or muon, and 
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For 771/ 3> 7770 we have the asymptotic behaviors 
^scalar — > 3/2 and £ 
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We begin with the electron system. As emphasized 
in ref. [7J, the (g — 2) e measurement is currently used 
to define the fine-structure constant a. The additional 
contribution to (g—2) e therefore acts as a shift of the fine- 
structure constant as Aa = 27rAa e . Comparing this cor- 
rection to measurements made in Rb and Cs atoms [HI [H] j 
the shift in a must not exceed 15 ppb, which constrains 
the coupling to electrons as [7J 

)\{m^ mi ) < 15 x 10" 9 . (7) 

For 7770 s» MeV this constraint translates to g e /e < 
2.3 x 10~ 4 and g e /e < 4.0 x 10 -4 for scalar and vec- 
tor mediators, respectively. The constraint is weakened 
for larger values of m$. 

More important for the purpose of a direct comparison 
with the Lamb shift in the /i-p system is the constraint 
coming from measurement of (<?— 2) M [TO]. At present, the 
theoretical prediction for a^f 1 seems to indicate a deficit 
of 302(88) x 10~ u compared with the experimental value 
Q exp jijj _ p or a given mass we can extract the values 
of g^ that bring the theoretical and experimental values 
of (g — 2) M into agreement. In fact, these values are in- 
sensitive to 7770 provided 7770 <C is satisfied, giving 
g^/e s» 1.6 x 10~ 3 for a vector and g^/e w 1.3 x 10~ 3 for 
a scalar. In Fig. [1] we plot the 2S 1/ / 2 -2P 3 / 2 energy shift, 
Eq. §, against the mass of the mediator, 7770, fixing the 
coupling to muons in this way. For the purpose of the 
plot we take g p = g M , but the result for other choices is 
easily obtained, as the energy shift is proportional to g p . 
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FIG. 1: The contribution to the energy shift in muonic hy- 
drogen, Eq. (pi, plotted against the mass of the mediator. 
In the central solid-blue curve we require the coupling to the 
muon, g M , to be such that the scalar contribution to (g — 2) M 
equals the theoretical deficit. In the upper/lower solid-blue 
curve the scalar contribution to (g — 2) M is determined to be 
±1 s.d. away from the theoretical deficit. The vector case is 
similarly given by the dashed-red curves. The coupling to the 
proton is fixed at g p = g^. The solid horizontal line is the dis- 
crepancy between the experimentally measured value of the 
energy and the theoretical prediction assuming the CODATA 
value [2] for the proton-radius, r p — 0.8768 fm. The dotted 
horizontal lines represent the ±1 s.d. uncertainty about this 
value. 



As the plot indicates, a new force with a mass s» MeV 
and a coupling to muons that explains the discrepancy in 
the muon anomalous magnetic moment can also give the 
required energy shift in muonic hydrogen to reconcile the 
proton-radius extracted from this system with the one 
extracted from hydrogen and electron-proton scattering. 
The mediator mass that gives the maximum energy shift 
is near an MeV because it is essentially determined by 
the Bohr radius of the fj,-p system, a^ 1 = 0.69 MeV. 
The choice g p s» g M favors ps MeV, but a different 
mass can be accommodated by increasing g p accordingly. 
For 7770 3> MeV (7770 <C MeV) the required coupling 
becomes large and scales as m 2 ^ (tt7^ 2 ). 

For 7770 MeV the coupling to muons necessary to 
explain both discrepancies turns out to be close to the 
muon mass divided by the electroweak scale, m^/v — 
4.3 x 10~ 4 . A Higgs-like coupling proportional to the 
mass would resolve any tension with the constraint from 
(g — 2) e on the electron coupling, Eq. However, it im- 
plies that the coupling to neutrons is comparable to the 
coupling to protons, g n ~ g pi which seems severely con- 
strained by neutron scattering experiments, as we discuss 
in the following section. 
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IV. CONSTRAINTS 

The most model-independent constraint on an MeV- 
scale interaction between muons and protons comes from 
measurements of the 3d 5 / 2 -2p3/2 transitions in 24 Mg 
and 28 Si [12] • A weighted-average between the re- 
sults from 24 Mg and 28 Si was used to obtain the limit, 
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Trifp Ri MeV this translates to a 95% CL limit g p g^/e < 
3.1 x 10~ 6 , assuming coupling only to protons. This mea- 
surement therefore allows the couplings necessary to ex- 
plain the muonic hydrogen and {g — 2)^ discrepancies. 
Mild tension does arise if one assumes that the new force 
carrier also couples to neutrons since the bound improves 
by a factor of 2. However, as we now discuss, such a 
coupling to the neutron is much better constrained by 
neutron scattering experiments. 

A new force carrier with an MeV mass that couples 
to neutrons produces sizable corrections to the scatter- 
ing cross-section of neutrons on heavy nuclei. This was 
first realized by Barbieri and Ericson [T2] who used the 
results of an old experiment [14] on the polarizability 
of the neutron to set a limit on any additional force car- 
rier that interferes with the strong interaction amplitude. 
They showed that such a force will contribute to the an- 
gular dependence of the differential cross-section in a dis- 
tinct manner as compared with the contribution from the 
strong interactions. Assuming that the new force couples 
equally to both protons and neutrons, the bound on the 
nucleon coupling is 

g n < 2 x 1(T 5 (m / MeV) 2 . (8) 

This is about an order of magnitude smaller than the nec- 
essary coupling to muons and protons discussed above. 
A more recent analysis arrived at a similar result |151 116] 
although the claimed precision of that experiment was 
later questioned (see Ref. [17] and references therein). 

The bound in Eq. ^ arises from an interference term 
between the strong amplitude and the amplitude due to 
the new force. As such, it may be susceptible to can- 
cellations involving other parts of the amplitude, and it 
also depends on the relative phase of the strong and new- 
physics amplitudes. Nevertheless, the constraint is suf- 
ficiently strong to disfavor g n w g p . The bound on the 
neutron coupling would have to be invalidated by more 
than an order of magnitude to allow for a simultaneous 
explanation of (g — 2) M and the muonic hydrogen results 
while requiring g n ss g pi a possibility which we therefore 
eschew in this letter. 

There are several other constraints on such a light bo- 
son, but they all involve further assumptions about its 
couplings to matter. Refs. [H1[T9] use 16 and 4 He atoms 
to search for + — > + + cf> transitions to constrain light 
bosons with Higgs-like couplings. These bounds are not 
very useful for ss MeV and in any case depend sen- 
sitively on the decay properties of <p. For example, the 
bounds do not apply if 4> decays promptly, or if <f> is too 
light to decay into an electron-positron pair. Ref. [50] 
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FIG. 2: The contribution to the 2S-2P energy splitting in 
muonic helium, Eq. (9j, plotted against the mediator's mass. 
The curve is normalized to SE^h = —0.31 meV, which is 
the discrepancy between the experimentally measured value 
and the theoretical prediction, assuming the CODATA value 
for the proton-radius, r p — 0.8768 fm. The dotted curves 
represent the ±1 s.d. uncertainty about this value. 



has searched for a light boson emitted in the scattering 
of electrons against the nucleus. It sets a strong bound 
on the coupling to the electron for masses in the range 
1.2 < ?Ti0 < 52 MeV. However, it too relies on the <f> 
lifetime being in a certain range. Ref. |21) sets a very 
strong limit on the decay ir + — > e + v e {4> — > e + e~) and 
excludes masses in the range 10 < < 110 MeV, but 
again relies on the coupling to electrons. Finally, we note 
that any coupling to neutrinos must also be strongly sup- 
pressed as such an interaction will strongly affect the well 
measured interactions of neutrinos with matter. 



V. CONTRIBUTIONS TO ENERGY SHIFTS IN 
OTHER MUONIC SYSTEMS 

We now move on to discuss predictions concerning 
other muonic systems such as muonic deuterium and he- 
lium, and true muonium. Neglecting any possible cou- 
pling to the neutrons, the energy shift in the 2S-2P tran- 
sition due to the new force is given by a simple gener- 
alization of Eq. |2]). Accounting for the change in the 
reduced mass and the atomic number, it can be written 
in terms of the contribution to the energy shift in muonic 
hydrogen, 



5E 



(jj,N) _ 




(9) 



where N stands for the different possible nuclei (deu- 
terium, helium, and etc.) and Z is the atomic number. 

The most straightforward prediction is that, more or 
less independent of the mediator's mass, the muonic deu- 
terium system, [i-D, should exhibit almost exactly the 
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FIG. 3: The contribution to the energy shift in true muonium 
as a function of the force carrier mass. Curves are as in Fig.[l] 

same energy shift in the 2S-2P transition as the /i-p sys- 
tem. Depending on m<*, it deviates from it by at most 
rs ±15% due to the change in the reduced mass, and we 
therefore predict 5E { £ D) = -0.3 ± 0.1 meV. 

Next we consider the muonic helium system. The 
2S 1 / 2 -2P 3 / 2 and 2Si/ 2 -2Pi/2 level splittings in 4 He) + 
muonic ion were measured in the 70's and reported in 
refs. [22l [23]. Unfortunately, these measurements were 
later criticized by other groups which cast doubt on the 
validity of the results. Experimentally [23], the 2Si/2- 
2P 3 / 2 resonance was not found at the frequency reported 
by ref. [22]. Theoretically [25], the lifetime of the 2S state 
is predicted to be too short at the pressure (40 bar) used 
in the experiment, a result which was later confirmed 
experimentally by ref. 20 . However, further measure- 
ments on muonic helium are planned in the near future 
at PSI pQ. In Fig. [2] we plot the prediction for the 2S-2P 
energy-shift in this system as a function of the mediator 
mass, ?n 0. 

A new force that couples to anions would also con- 
tribute to the energy levels of true muonium 
bound state), a system that is yet to be observed. In an 



ineresting recent proposal, the authors of Ref. [27] have 
discussed two possible production mechanisms, e + e~ — > 
and e + e~ — > 7/j + /i~, as particularly promising 
channels. If in addition to production, spectroscopic 
studies of true muonium can be achieved, then it may 
be possible to observe the small energy-shift due to the 
force carrier exchange. In Fig. [3] we plot the predicted 
energy shift in the 2S-2P transition when the muon cou- 
pling, <7 M , is chosen so as to fix the discrepancy in (g— 2)^. 
This energy shift is much smaller than the expected level- 
splitting due to vacuum polarization and will undoubt- 
edly be very difficult to observe. Nevertheless, it serves 
as an unambiguous prediction that may be tested in the 
future. 

VI. CONCLUSIONS 

A new force carrier with a mass of ~ MeV that cou- 
ples to both protons and muons with g^ :P ss 4 x 10 -4 can 
explain the discrepancies observed in both (g — 2)^ and 
the muonic hydrogen 2S!/ 2 -2P3/2 energy splitting. How- 
ever, the coupling of such a new force to cither neutrons 
or electrons is constrained by several past measurements. 
Nevertheless, none of these constraints exclude this pos- 
sibility in a model-independent way. Assuming only cou- 
plings to muons and protons, this possibility lends itself 
to concrete predictions for the expected energy shifts in 
muonic deuterium and helium as well as true muonium. 

Note added: While this paper was being completed 
ref. [28] appeared, which also considers new-physics con- 
tributions to the Lamb shift in muonic hydrogen. 
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